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We study the production of photons through the non-equilibrium relaxation of a disoriented chiral 
condensate formed in the expanding hot central region in ultra-relativistic heavy-ion collisions. It 
is found that the expansion smoothes out the resonances in the process of parametric amplification 
such that the non-equilibrium photons are dominant to the thermal photons over the range 0.2-2 
GeV. We propose that to search for non-equilibrium photons in the direct photon measurements of 
heavy-ion collisions can be a potential test of the formation of disoriented chiral condensates. 

PACS numbers: 25.75-q, 11.30.Rd, 11.30.Qc, 12.38.Mh 

In relativistic heavy ion collisions, the highly Lorentz contracted nuclei essentially pass through each other, leaving 
behind a hot plasma in the central rapidity region with large energy density corresponding to temperature above 
200 MeV where the chiral symmetry is restored. This plasma then cools down via rapid hydrodynamic expansion 
through the chiral phase transition during which the long-wavelength fluctuations become unstable and grow due to 
the spinoidal instabilities [Q, ||, || . The growth of these unstable modes results in the formation of disoriented chiral 
condensates (DCCs) (see also Refs. j| ||). The DCCs are the correlated regions of space-time where the chiral order 
parameter of QCD is chirally rotated from its usual orientation in isospin space. Subsequent relaxation of the DCCs 
to the true QCD vacuum is expected to radiate copious soft pions which could be a potential experimental signature 
of the chiral phase transition observable in the ultra-relativistic heavy ion collisions at RHIC and LHC although such 
a signature has been looked for and not seen in the lower energy heavy ion collisions at the CERN SPS [|. However, 
since these emitted pions will undergo the strong final interaction, the signal may be severely masked and become 
indistinguishable from the background. It then becomes important to study other possible signatures of DCCs that 
would be less affected by the final state interaction. Electromagnetic probe such as photon and lepton with longer 
mean free path in the medium serves as a good candidate and can reveal more detailed non-equilibrium information 
on the DCCs with minimal distortion jj], |§|, Hj. 

Minakawa and Muller jlCj have suggested that the presence of strong electromagnetic fields in relativistic heavy 
ion collisions induces a quasi-instantaneous "kick" to the field configuration along the 7r° direction such that it is 
plausible that the chiral order parameter in the DCC domains, if formed, will acquire a component in the direction 
of the neutral pion. The production of photons through the non-equilibrium relaxation of a DCC within which 
the chiral order parameter initially has a non-vanishing expectation value along the ir° direction and subsequently 
oscillates around the minimum of the effective potential has been considered. In Ref. ||, Boyanovsky et al. have 
extensively studied the photon production from the low energy coupling of the neutral pion to photon via the Ua(1) 
anomalous vertex. They have found that for large initial amplitudes of the 7r° field photon production is enhanced by 
parametric amplification. These processes are non-perturbative with a large contribution during the non-equilibrium 
stages of the evolution and result in a distinct distribution of the produced photons. Later, the authors in Ref. pl| ] 
have taken into account another dominant contribution that also involves the dynamics of ir° due to the decay of 
the vector meson through the electromagnetic vertex. Although the corresponding dimcnsionless effective coupling 
involving the vector meson is quite perturbatively small, for the large amplitude oscillations of the 7r° mean field, 
the contribution to the photon production is of the same order of magnitude as the anomalous interaction. However, 
they have ignored the hydrodynamical expansion and adopted the simple "quench" phase transition from an initial 
thermodynamic equilibrium state at a temperature (Tj) higher than the critical temperature (T c ) for the chiral phase 
transition cooled instantaneously to zero temperature, which has been widely used in the study of non-equilibrium 
phenomena of DCCs §, @, |, §. 
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In this paper, we will consider the effect of the hydrodynamical expansion of the plasma to the production of 
photons from the non-equilibrium relaxation of a DCC. Although the "quench" scenario can capture the qualitative 
features of this non-equilibrium problem, we will see that the dynamics of parametric amplification as well as spinoidal 
instabilities is considerably modified by the plasma expansion. In fact, our study of the photon production based 
upon this strongly out-of-equilibrium scenario bears analogy to the recent studies in the context of the off-shell 
photon production from the out-of-equilibrium expanding quark-gluon plasma with a finite life time [p^| . To make 
a comparison, we also compute the thermal equilibrium photon emission from the hot quark-gluon plasma as well 
as the hadronic matter following an adiabatic quark-gluon phase transition by convoluting the equilibrium photon 
production rates at finite temperature with the expansion dynamics. 

In ultra-relativistic heavy-ion collisions, it is expected that the rapidity density of the particles produced in the hot 
central region has a plateau. This implies an approximate Lorentz boost invariance along the longitudinal direction in 
the evolution of the hot plasma in the central region jL3|. However, at late times following the heavy nuclei collisions, 
a transverse flow can be generated due to the multi-scattering between the produced particles, as such the expansion 
becomes three dimensional |Q. Here we will simply assume that the hydrodynamical flow is spherically symmetric, 
and that the boost is along the radial direction. We will see that this assumption greatly simplifies the treatment of 
the photon field in the expanding space-time. We consider this as the first simplest attempt to tackle the problem. 

The natural coordinates for spherical boost invariant hydrodynamical flow are the proper time r and the space-time 
rapidity 77 defined as 

r = (i 2 -r 2 )i ^=1 ]n (111), (1) 

where (t, r) are the coordinates in the laboratory, 

t = TCOsh?7, r^rsinhr;. (2) 

The ranges of these coordinates are set to be < r < 00 and < 77 < 00, restricted to the forward light cone. In 
terms of spherical coordinates, the Minkowski line element is given by 

ds 2 = df-dr 2 (3) 
= d^-T^d^ + sinh^rf^+sinh^sin 2 ^^ 2 ). (4) 

Eq. (|4|) is the Robertson- Walker metric for an open expanding universe with the scale factor R(t) = r Jlf|. For 
77 < 1, the open metric (|J) can be approximated by the flat Robertson- Walker metric, 

ds 2 = dr 2 - T 2 {dr] 2 + j] 2 d6 2 + j] 2 sin 2 d(/) 2 ), 

= dT 2 -r 2 drf. (5) 

This flat metric will serve as the background comoving frame under which to study non-equilibrium photon production 
from DCC domains. This is a good first-order approximation since the boost invariance is in anyway applied only 
for small values of 77. Furthermore, 77 < 1 corresponds to the space-time with r/t < 0.76 which has covered a major 
portion of the forward light cone. 

The relevant phenomenological effective action in a general expanding space-time is given by 

S = J d A Xyfg (La- + L A + L n o A ) , (6) 

where 



1 M 

L a = ~ff^5 M *-ft,$+--^*-$-A(*-$) +ha, (7) 



La = -\g^g^F a& F^ (8) 

L^oa = -l-£L^ e °^ F f + 1 ffv dx7: d F - F (9) 

where F M „ = d^A v — d v A^ is the photon field, and $ = (a, tt°,tt) is an 0(N + 1) vector of scalar fields with 
7? = (n 1 , 7r 2 , -n- 7 ^ -1 ) representing the N — 1 pions. Note that the signature is ( — h ++), and we have added 
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to each e al3tll/ because it is a tensor density of weight —1 |T^]. The phenomenological parameters in the effective 
Lagrangian above can be determined by the low-energy pion physics as follows: 



m a fa 600 MeV, U 
h w (120 MeV) 3 



93 MeV, A Pa 4.5, T c Pa 200 MeV, 
m v « 782 MeV, \ v ps 0.36, 



(10) 



where V is identified as the ui meson, and the coupling Ay is obtained from the u> —> tt°j decay width [ fl6| . This 
effective action has been obtained by two of us jllj] in the study of out-of-equilibrium photon production from DCC 
domains by taking the rapid "quench" phase transition scenario. Here we will consider the similar non-equilibrium 
phenomena taking account of the hydrodynamical expansion. It must be noticed that this is an effective field theory 
with an ultraviolet momentum cutoff of the order of A sa 1 GeV. The effective vertices obtained are in terms of 
the perturbative theory without involving in-medium modifications. In-medium effects will enter only through the 
non-equilibrium Green's functions for the meson fields as we will see below. To obtain these effective vertices including 
the strongly out of equilibrium effects, one should integrate out the quark fields and the vector meson in the context 
of the fully non-equilibrium formalism. It deserves to tackle in the near future. 

It is well known that the minimal coupling of photons to the metric background is conformally invariant. Thus, in 
the conformally flat metric (0), it is convenient to work with the conformal time defined by 



du = 



(It 



(11) 



where Tj is the initial proper time after which we expect that the quark-gluon plasma is formed and thermalized. 
Integration gives 



n In | — 

Ti 



(12) 



Rewrite the flat metric 



as 



ds 2 = — g^ v dx^dx v — a 2 (u)(du 2 — dx 2 ), 



(13) 



where a(u) — t/tj = e u ^ Ti and dx = Tidff. Writing explicitly in terms of the comoving coordinates and defining 
<f> = $ a /a, the action (0) becomes 



where 



Co 



S = J du d 3 x C = J du d 3 x (C a + Ca + C-k°a) > 



1 



.'I '/ r uf3 r [Li>i 



e 2 tt° 



1 9 , ,o 1 d 2 a 
2 aM ° + a^ 



(*a'*a)' 



a ha a , 



32ir 2 aU 



a 2 8m 2 my 



(14) 

(15) 
(16) 
(17) 



where $ a = (er a , 7r°, ff a ) and rf v is the Minkowski metric. In terms of the conformal time, the effective action now 
has analogy with the effective action in Minkowski spacetime with the time dependent mass term and interactions. 

The dynamics of the non-equilibrium expectation values as well as correlation functions of quantum fields can 
be obtained by implementing the Schwinger-Keldysh closed-time-path formulation of non-equilibrium quantum field 
theory. This formulation is used to describe the evolution of an initially prepared density matrix, and requires a path 
integral defined along a closed time contour. This techniques have successfully been employed elsewhere within many 
different contexts and we refer the readers to the literature for details p7| . Here we first shift a a and 7r° by their 
expectation values with respect to an initial non-equilibrium states: 



cr Q (x, u) = (f> a (u) + Xa (x, u) , 

u ) = Ca(u) + l/j a (x, u) , 



(<Ta(x, u)) = 4>a(u), 
(^aO^")) = Ca(u), 



(18) 
(19) 
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with the tadpole conditions: 



(Xa(x,u)) = 0, (tp a (x,u)) = 0, (Tt a (x,u)) = 0. 



(20) 



One can impose these tadpole conditions to all orders in the corresponding expansion to derive the non-equilibrium 
evolution equations of field expectation values. In addition, the large-iV approximation will be implemented below 
(Eqs. (2.7)-(2.9) of Ref. 0) to obtain the factorized Lagrangian which provides a non-perturbative framework to self- 
consistently incorporate quantum fluctuation effects from the strong a — ir interactions. As for the issue of the gauge 
choice, we adopt here the Coulomb gauge which involves only the physical degrees of freedom without any redundant 
fields to study non-equilibrium photon production from the oscillating mean fields within the DCC domains in a 
spherical expansion. 

Using the tadpole conditions (po|), the large- TV equations of motion for the mean fields involving the perturbatively 
electromagnetic corrections can be obtained as follows: 



du 2 

du 2 
+- 



+ m 2 a (u) + 4\<j>l(u) + 4\£(u) + 4A(t?2)( u ) 
+ ml(u) + 4A#(u) + 4\£(u) + 4A(t?2)(u) 



a 3 h = 0, 



327r 2 / T a 



(FF) 



e 2 X 2 v 1 d fCa 



a 2 du 



Am^rriy a 3 du \ a 



Vasid-yF^Fap) = 0. 



We then decompose the fields 7? a and At into their Fourier mode functions U ak *(u) and V Xk *(u) respectively, 

d 3 k 



n a (x, u) 
A T (x, u) 



2(2tt) 



d 3 k 



d 3 k 



-.A T (k,u) 



(21) 



(22) 



{ [b + kV lZ (u)e +j : + b_^V 2S (u)e_ s 



h.c 



(23) 



where a^ and 6 ± £ are destruction operators, and e* ± g are circular polarization unit vectors. The frequencies g i 
and u> A t . can be determined from the initial states and will be specified below. The mode equations for the photons 
as well as the pions can be read off from the quadratic part of the effective factorized Lagrangian as 



du 2 

' 

du 2 

' £_ 

du 2 



- ml(u) + i\<t> 2 a (u) + 4\£(u) + 4X(7f 2 a )(u) 

-C^j v lh = o, 

<k I V 2k = 0, 



Uak = 0, 



e 2 \ 2 v 


c 2 


k 2 - 


e 2 


9 2 


a 2 


2w 2 U 


e 2 \ 2 v 


c 2 " 


k 2 + 


e 2 


mlm 2 v 


a 2 


2tt 2 U 



(24) 
(25) 

(26) 



where m 2 (u) = —(1/2) a 2 M 2 — a/a, ((u) = ( a (u)/a, k — \k\, and the dot means d/du. The expectation values with 
respective to the initial states are given by 



e a ^(F af3 F^)(u) 
^ 0a e am ?l*8(F^F a p)(u) 



(N-l) 

Aa(u) 



Aa(u) 



d 3 k 



2(2^) 3 w ffafe 



[\U ak (u)\ 2 ] coth 



22} 



5 ^_ (4t) | nWlW |.-|« lW | ! 



Aa(u) 



d 3 k 



2(27T) 3 UUfc 



(4k 2 )[\V lk (u)\ 2 + \V 2k (u)\ 2 ], 
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Aa(u) 



,r ' k -m^[\v lk (u^ + \v 2k (u)\ 2 ], 



where (7?^) (it) is self-consistently determined by Eq. 



(^)(0) = (iV-l) 



2(2Tt) 3 LUAk v ' du 
4J). In particular, when u = 0, it becomes the gap equation 
A d 3 k 



(27) 



2(2^) 3 ^ afe)i 



coth 



2Ti 



(28) 



Notice that the photonic medium effects have been ignored in the above expressions of the photon field expectation 
values due to the fact that the photons interact electromagnetically so that their mean free paths are expected to be 
longer than the estimated size of the quark-gluon plasma fireball, and the produced photons will escape freely toward 
the particle detector having no enough time to build up their population in the plasma. In this conformal frame, 
the momentum cutoff we choose depends linearly on a so as to keep the physical momentum cutoff A fixed in the 
laboratory frame. As a result, the initial conditions for the mode functions become more subtle here. Here we choose 



U a k<A{0) = 1, ^afe<A(0) = -«Wff a fc,i, 

c4 oM = fc 2 + m 2 (0) + 4A^(0) + 4AC 2 (0) + 4A(t? 2 ,)(0); 
Vi,2fc<A(0) = 1, Vi,2fc<A(0) = -iu>Ak,i, UAk,i = k; 

U a k>k{uk) = 1, U ak> A(uk) = -i&i? a k,i, 

<4 B k,< = fc2 + m aK) + ^l(uk) + 4AC 2 K) + 4A(7?2)( Mfc ); 
Vi,2fc>A(Wfc) = 1) Vl,2fe>A(wfc) = -MAk,i, U A k,i = k. 



(29) 



(30) 



Therefore, for the momentum lying below the momentum cutoff A at the initial conformal time, the initial condition 
for the mode function can be set at u — as shown in Eq. (^). However, for the momentum above the cutoff 
momentum A at u = 0, it will become dynamical when k < Aa(u) after u = u k = rdn(fc/A), and the initial condition 
for this mode function must be set at that conformal time as in Eq. ( |30| ) . The photon spectral number density at 
time u is given by the expectation value of the number operator for the asymptotic photons H , 



dN 



dN 4 



dN_ 



d 3 xd 3 k d 3 xd 3 k d 3 xd 3 k 
= (N fc (u)> 
1 

2k 
1 

4P 



A T (k, u) ■ A T (-k, u) + k 2 A T {k, u) ■ A T (-k, u) 



- 1 



\V lk (u)\ 2 + k 2 \V lk (u)\ 2 
and the invariant photon spectral production rate is given by 
kdN 



1 

sLIO 



\V 2k (u)\ 2 + k 2 \V 2k (u)\' 



e 2 \ 2 v 



dud 3 xd 3 k 



\V 2k \' 



C- 



2ir 2 f 7T du 



\Vu 



(31) 



(32) 



We now need to relate this spectral production rate to the invariant production rate measured in the metric in 
order to obtain the photon spectrum in the laboratory frame. From the coordinate transformations (|^), we find that 

q = p cosh 77 + p r sillily, 
q v = pr sinh r) + p r T cosh 77, 
qe = Pe, 

q<i> = p<t>, 



(33) 



and the zero- mass condition of the photon, i.e., g^ v p^p v = 0, 

1 



q = M = I 
r r 



1 



2 

sinh 7] 



\p\ 



1 



Pl + —Pe 



1 



sinh ry sin 9 



(34) 
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where — (q, q) is the photon four-momentum in the metric (Q) and p^ = (p,p) is that measured in the laboratory 
frame. Because of the spherical symmetry of the problem, for convenience, we can choose p to be along the z-axis, 
which is also the polar axis. This gives p r = pcos9, and then from Eq. (p3|), 



q = p(cosh?7 + cos sinh??). (35) 



Hence, from the momentum transformation laws (p3[), one can change the momentum variables from q^ into p^ in 
the invariant production rate making use of the fact that the spectral particle number density dN is invariant under 
such transformations. We then find that 

dN = gdN = pdN 

dT dq^dqedq^drd-qdOdcj) d 3 p r 3 sinh 2 r] sin 6 dTdrjdOd^' 

where dT stands for the invariant phase space element in the comoving frame. Therefore, the photon spectrum 
measured in the laboratory at the final time Tf can be obtained by integrating the invariant spectral production rate 
over the space-time history of the hydrodynamical evolution in the comoving frame as 

„, iA l Tf dr [ Vm ~ drj f d9 [ d(j> t 3 sinh 2 rj sin 0^- (t, q) , (37) 



d 3 P J Ti Jo Jo Jo dT 

where q — p(cosh?7 + cos#sinh?7). Now, in the small rj approximation, one can link the invariant rate in the comoving 
frame to that of the conformal frame as given in Eq. (|32]), where k = qr/ri and k = q/ri, in the following way: 

dN qdN /r 2 \2 kdN 

- ' ! (38) 



dT drd 3 ffd 3 q \tJ dud 3 xd 3 k 
where u — rdn(r/Ti). Then, we can approximate 

p dN n , /•*»— , f w jn f 2 * J± 2 . - .kdN 



d 3 p 



/•Tf fiBa /* f ■ k dN 

/ dr drj d0 d(j> r 2 Ti sinh 77 sin # -(u,k) 

jTi Jo Jo Jo dudxdk 

.^max t rr: r -1/7 

1 / KTi p T 



2 7rr 2 PTf l'TpB r > mac /T i 

dr / dk 



P J Ti J Tpe 7 ?max I 



cosh(?7 m ax) . , 

2 \p T 



k dN . , . 

-,{u,k) . (39) 



dud 3 xd 3 k 



The final expression is obtained after changing the integration variable from 9 to k as well as carrying out the <f>- and 
■q- integrals. 

In contrast with the non-equilibrium photon production from the DCC domains, the equilibrium photon emission 
arising from an expanding quark-gluon plasma (QGP) as well as hadronic matter can be obtained by convoluting 
the equilibrium photon production rates at finite temperature with the expansion dynamics. We currently believe 
that after the ultra-relativistic heavy collisions, in the central rapidity regime, parton-parton multi-scatterings cause 
the produced quarks and gluons to reach local thermal equilibrium at the time scale of Tj, and then the so-called 
QGP with high energy/entropy density is expected to be formed. After the formation of a QGP, the dynamics of 
the thermal plasma is dominated by the hydrodynamical expansion. The hydrodynamical approach is appropriate 
since in this stage quarks and gluons interact in such a dense regime that their mean free paths are much smaller 
than typical wavelengths of various collective phenomena. Subsequently, the plasma cools adiabatically according to 
the equation of state of a hot QGP down to the quark-hadron phase transition when the QGP starts to hadronize. 
Assuming that the phase transition is a first order one, quarks, gluons, and hadrons will coexist in a mixed phase in 
which the released latent heat keeps the temperature of the plasma unchanged at the critical temperature even though 
the plasma continues to expand. This mixed phase will exist until all the matter has converted into the hadrons. The 
system in the hadronic phase will continuously expand with the dynamics governed by the equation of state of a hot 
hadron gas. As a result, the temperature continuously drops till the so-called freeze-out temperature is reached where 
the mean free path for all hadrons is of the order of the size of the plasma, and after that all hadrons cease further 
interactions. We would like to emphasize at this stage that according to such a hydrodynamical picture, the strongly 
out of equilibrium scenario of either supercooling or superheating has been ignored in the sense that the entropy 
of the system is assumed to be locally conserved for the whole evolving history of the system. Therefore, photon 
emission under this adiabatic expansion is in striking contrast with that from the DCC domains associated with a 
strongly out of equilibrium scenario. In the Bjorken's hydrodynamical model, the QGP as well as hadronic matter 
can be described as a fluid in local thermal equilibrium. The energy density, pressure, entropy etc. will depend on the 
proper time only, and will not depend on the space-time rapidity under the assumption of a boost invariance in the 
central rapidity region. The dynamics of the fluid is determined by the conservation laws of the energy-momentum, 
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baryon number, and entropy. In addition, an equation of state of the fluid in the various phases can serve as an input 
to close this set of dynamical equations to be solved self-consistently. Here, to be consistent with the calculation 
of the non-equilibrium photon production from DCC domains, we also assume that the underlying hydrodynamical 
expansion is of a spherical expansion with the coordinates of the proper time r and the space-time rapidity r\ that 
results from free-streaming particles with a collective velocity v — r/t along the radial direction, where the fluid 
rapidity is equal to the space-time rapidity. In the central rapidity region which is baryon free, the conservation law 
of entropy for such an adiabatic expansion reads 

s(T)t 3 = s(T t )rl (40) 

This conservation law will play an important role in determining the cooling law as well as the duration time scales 
of various phases of the evolving system from the initial QGP phase, to the quark-hadron mixed phase, and to the 
final hadronic phase. 

In the QGP phase, the entropy density is dominated by the relativistic massless gas composed of quarks and gluons 
given by 

sq = 2 -^9gT\ (41) 

where the degeneracy go = (2 • 8 + 2 • 2 • 3 • 2 • |) for the two-flavor quarks and the eight SU(3) gluons. On the other 
hand, the entropy density of hadronic matter is taken as 

sh = ^W)T 3 , (42) 

where gn{T) is an effective degeneracy for a hadron gas. We introduce a temperature dependence in this effective 
degeneracy to account for the medium effects from the coupled hadron gas in the temperature regime of our interest. 
The variation of the effective degeneracy can be parametrized as gn(T) ~ 13(T/fm -1 )' 5 where the generally acceptable 
value of S is about S ~ 3.4 of slight model-dependence (here we adopt the Walecka model) |lq ]. We are now in a 
position to determine the time scale in the QGP phase which starts from the initial formation time t\ until tq when 
the temperature drops to the critical temperature T c of the quark-hadron phase transition and a phase transition 
to hadronic matter starts. Following the conservation law of entropy together with an expression of sq in Eq. (|4l|), 
the cooling law in the QGP phase can be obtained as T = (rj/r)Tj, and thus tq = (Tj/T c )rj. However, during 
the quark-hadron mixed phase, the relative proportion of QGP to hadronic matter must vary with time. Since the 
temperature of the quark-hadron mixed state remains constant at T = T c even though the plasma continuously 
expands, the entropy density of the mixed state is equal to the sum of the entropy densities of QGP and hadronic 
matter at T = T c , Sq and s H , respectively weighted by their fractions, /q(t) and /h(t): 

Smix(r) = Sq{t)s c q + f H (T)s c H , (43) 

where /q(t) + /h(t) — 1- Initially at r = tq, the phase consists entirely of QGP, i.e. /q(tq) = 1 , and at the end, 
entirely of hadronic matter at th, i.e. fH{ T H) = 1- Thus, the life time of the mixed state is t# — tq. /q(t), /h(t), 
and th are to be determined from the entropy conservation law (^0|) which presumably also holds in the mixed state. 
Then, from Eqs. (§||) and (|(j|), we have 

(/q(t> c q + [1 - !q(t)] Sh ) t 3 = ^r 3 . (44) 

Hence, we obtain 

\3 



fQ(r) = SQ[ c T ' - SH , /h(t) = 1-/q(t). (45) 



'TQ - 

S Q ~ °H 

>V3, 



Therefore, the mixed phase ends at th = ( s q/ s h) T Q determined from /q(t^) = above. As for the hadronic 
phase, substituting the associated entropy density of a hot hadron gas ( pl2| ) into the entropy conservation law (|4(]) can 
lead to the cooling law in this phase given by T = T , c (th/t) 3 /' 5+3 - ) from which the freeze-out time can be obtained as 
T f = T„(T c /T f )V+W. 

After having obtained all relevant duration time scales in various phases from the QGP phase to the hadronic 
phase, we now turn into the discussions of the photon production processes from the QGP and the hadronic phase. 
The calculation of the thermal photon production rates in the high energy regime (E >> T) from the QGP at 
one-loop order in the hot thermal loop (HTL) approximation has been obtained by Kapusta et al. to account for 
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the real photon production from the annihilation of a quark-antiquark pair into a photon and a gluon (qq — > gj) as 
well as the absorption of a gluon by a quark (antiquark) emitting a photon (q(q)g — * 9(9)7), similar to Compton 
scattering in QED JTsj] . However, the recent study by Aurenche et al. has shown that the two-loop contributions 
to the photon production under the HTL approximation from quark (antiquark) bremsstrahlung ( 99(g) — > 99(3)7), 
and quark-antiquark annihilation with scattering ( qqq(g) — * 9(5)7) are of the same order of magnitude as that in 
one-loop order fl20fl . Their contributions to the photon production rates prove to be dominant at high photon energies. 
Later, we will adopt this two-loop result (Eqs. (1,2,3) of Ref. |2l|) to calculate the equilibrium photon emission during 
the phases that involve the QGP. As for the photon production processes from the hadronic phase, it is known that 
the dominant processes in the energy regime of our interest are those of the reactions tttt — > p7, np — > as well as 
the decays ui — ► ttJjP — * ttttj p2| . The reaction for producing photons that involves the intermediary axial vector 
meson ai is also important, for example, 7rp — > ai — > 777 |23fl . All of these hadronic processes will be taken into 
account in the calculation of the equilibrium photon emission from the thermal hadrons. A point to be stressed is 
that we consider the photon production processes from the low-mass hadrons only while the higher mass resonances 
would be Boltzmann suppressed. We then use the invariant photon production rates either from the QGP or from the 
hadrons combining with the relevant cooling laws as well as their expansion dynamics to numerically obtain the total 
equilibrium photon emission for the whole evolution of the system. Notice that the equilibrium photon production 
invariant rates obtained in the literature fll^ , Q in the rest frame of the emitting matter can be expressed in terms 
of the comoving frame of a thermal bath with temperature being a function of the proper time r under the adiabatic 
expansion. We can denote those invariant rates from the QGP and the hadron gas in the comoving frame respectively 
by dN/dT \q and dN/dT \h that both depend on the photon four-momentum g M as well as the temperature T(t) 
with r dependence following their respective cooling laws. Finally, the equilibrium photon emission starting from the 
formation time Ti down to the freeze-out time 77 in the laboratory frame can be obtained by a Lorentz boost from 
the comoving frame to the center of mass frame. Thus it becomes with Eq. ( [37| ) 

rl AT r T Q PVmELX pit p2tt rl AT 

^ = J dr J d V J dO J d4> r 3 sinhSsin^ \ Q (q,T = 
r™ . . r „ f 2 * . , , . , / , , diV dN 



dr ^ d V ^ dO ^ d<f> ^sinhVinfl ^/q(t)- \ q {q,T c ) + f H (T)— \h {q,T c ) 
dr J dr] J d6 J d<f> r 3 sinh 2 ^sin6»— \ H [q,T = T c (— J 



(46) 



where again q = p(coshr/ + cos 9 sinh if) in Eq. (35). 

We now perform the numerical study. We first solve the coupled differential equations Eq. (^l|), Eqs. ( p4[ - p6] ), and 
the field expectation values Eq. ( p7[ ) self-consistently with the initial conditions £(0) = 0.3 fm~ and £(0) = 0(0) — 
</>(0) = in the conformal time frame. The initial temperature Ti at the QGP formation time Ti = 1.0 fm _1 that 
corresponds to the initial conformal time u — is chosen to be Ti = 1.0 fen . We set N = 3 for 3 pions. In Figs. 1(a) 
and 1(b), we show that C( u ) evolves with strong damping, while the amplitude of the field derivative C( u ) grows. The 
strong damping is due to the back-reaction effects from the creation of excitations of the quantum fluctuations via 
parametric amplification. This translates into photon production as we will discuss later. In fact, the growing £(w) 
makes photon production become more effective. The spectra of emitted photons in the conformal time frame are 
depicted for u — 1, 1.5, and 2 fm^ 1 respectively in Fig. 1(c). They feature prominent peaks. Notice that the peaks 
are located at k ~ w^/2 and k ~ u>£, where is the time-averaged oscillation frequency of C which can be read off 
from Fig. 1(b). Since uj^ increases with time, both of the peaks move towards the higher momenta as we can see from 
Fig. 1(c). Therefore, it is evident that the photon production is due to parametric amplification. Figs. 2(a), 2(b), 
and 2(c) are shown for a large initial amplitude £(0) = 1.0 fm _1 and C(0) = 0(0) = 0(0) = 0. This is clearly the 
non-perturbative effects driven by the parametric amplification which exhibits the features of unstable bands leading 
to the growth of the fluctuation modes. The growth of the fluctuations turns into the profuse photon production in 
the modes within the unstable bands. It is not surprising that with a large initial field expectation value, the peaks 
in the photon spectrum (Fig. 2(c)) are about an order of magnitude larger than the previous case. 

The photon spectrum in the laboratory frame emitted from the non-equilibrium DCC domains can be obtained 
by integrating the rates ([52|) over the spacetime history (|3^). To make the comparison, we calculate the thermal 
photon emission from the adiabatic quark-hadron phase transition in Eq. (Q) with the photon production rates for 
the processes discussed above. To do so, as for the non-equilibrium DCC case, we assume that the QGP is formed 
with temperature Ti — 1.0 fm -1 at the formation time Tj = 1.0 fen. In addition, the cooling law in the QGP phase, i.e. 
t = {T/Ti)Ti, gives the time scale tq = 1.1 fm when the QGP begins hadronizing and the subsequent quark-hadron 
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FIG. 1: (a) Temporal evolution of the 7r° mean field ({u). (b) Evolution of the derivative ((u). (c) Photon spectral number 
density dN /cPxdPk for £(0) = 0.3 fm -1 at different times, where u is in units of fm. 
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FIG. 2: As in Fig. 1 for C(0) = 1.0 fm" 1 . 
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FIG. 3: Photon energy spectrum measured in laboratory. The solid line denotes the initial mean field C(0) = 1.0 fin" 1 , while 
the dashed line denotes ("(0) = 0.3 fm -1 . The dotted line is the thermal photons emitted from quark-gluon phase transition. 



mixed state is formed at the critical temperature assumed to be T c = 0.9 fm -1 . One can also estimate the time when 
the mixed state ends, t h = (sq/s^) 1 / 3 = 1.77 fm, where the entropy densities are evaluated at T = T c . In the final 
stage, the system turns into the hadronic phase starting at th till the freeze-out time t/ = 5.6 fm estimated from 
the cooling law of the hadronic phase, ry = tjj (T c /TfY s+3 ^ 3 , where 5 ~ 3.4 jTgj and the freeze-out temperature is 
chosen as Tf — 0.5 fm~ . 

Figs. 3(a), 3(b), and 3(c) display the photon emission in the laboratory frame from the non-equilibrium DCC 
domains (solid line for £(0) = 1.0 fm -1 and dashed line for £(0) = 0.3 fm -1 ) and the adiabatic evolution of the QGP 
phase down to the hadronic phase (dotted line) for the spacetime rapidity rj — 0.65, 0.8, and 1.0 respectively. In all 
three cases, our results indicate that the non-equilibrium photon emission from DCC domains for an initial value of 
the 7T° component of an order parameter about 1.0 fm -1 that undergoes large amplitude oscillations with respect to 
the minimum value of the effective action is an order of magnitude larger than the photon emission from the QGP to 
the hadronic phase through an adiabatic first order phase transition. On the other hand, for an initial value of the 
7T° below 0.3 fm -1 , due to small amplitude oscillations, the non-equilibrium DCC photon emission is not so effective 
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as compared to the thermal photon emission. Therefore, we believe that if DCC domains are indeed formed, the 
nonlinear dynamics of the ir° oscillations can create an enormous amount of non-equilibrium photons from the DCC 
domains. These photon emission enhancement effects could potentially provide a distinct experimental signature. In 
addition, the emitted non-equilibrium photon spectra from DCC domains in a hydrodynamical expansion shown in 
Fig. 3 reveal the dramatically different features as compared with the results in Refs. M O, where the authors have 
considered the similar non-equilibrium phenomena in the temperature quenching scenario without convoluting the 
photon production rate with the plasma expansion, and shown that prominent peaks appear in the photon spectrum 
locating at resonant momenta of the parametric amplification determined by the frequency of the ir° oscillations. 
Since the measured emitted photons are chosen with the momentum along the z direction and the plasma expansion 
is along the radial direction, the Doppler effect would cause the emitted photons cither blue-shifted if the photon 
comes from the place where the direction of its momentum is along with that of the expansion at that place, or 
red-shifted if the direction of its momentum is opposite to that of the expansion. As a result, as we can see in Fig. 3, 
the Doppler effect leads to smoothing out the resonant peaks found in Refs. |9|, [fl|. In addition, one can also estimate 
the corresponding boost factor arising from the background expansion effect. From Eq. ( |35| ) and k = qr/n, one 
can relate the momentum k in the comformal frame to the momentum p in the laboratory frame by p ~ e^fc, where 
rj = 0.65, 0.8 and 1.0 in Fig. 3. This background expansion may boost the photon energy, say below the ultraviolet 
momentum cutoff of 1 GeV produced from the low-energy effective interactions in Eqs. (0^), to the energy as high 
as the order of 2 GeV. 

In conclusion, we have studied the production of photons through the non-equilibrium relaxation of a disoriented 
chiral condensate within which the chiral order parameter initially has a non-vanishing expectation value along the 
7T° direction, taking into account the hydrodynamical expansion of the central rapidity region. The production of 
non-equilibrium photons driven by the oscillation of the 7r° field due to parametric amplification is modified by the 
expansion which smoothes out the resonant photon production in the "quench" scenario. However, the resulting 
photon spectrum exceeds that for thermal photons from quark-gluon plasma and hadronic matter for photon energies 
around 0.2 — 2 GeV. These non-thermal photons can be a potential test of the formation of disoriented chiral 
condensates in relativistic heavy-ion-collision experiments. 

We would like to thank D. Boyanovsky, Hector de Vega, and S.-Y. Wang for their useful discussions. The work 
of Y.-Y. Charng, K.-W. Ng, and D.-S. Lee were supported in part by the National Science Council, ROC under the 
grants NSC90-2811-M-001-071, NSC90-2112-M-001-028, and NSC90-2112-M-259-011 respectively. 
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